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Abstract 

Visual and quantitative information were obtained 
on heat transfer and flow in a branched-duct test section 
that had several significant features of an internal cool- 
ing passage of a turbine blade. The objective of this 
study was to generate a set of experimental data that 
could be used to validate computer codes for internal 
cooling systems. 

Surface heat transfer coefficients and entrance flow 
conditions were measured at entrance Reynolds numbers 
of 45,000, 335,000, and 726,000. The heat transfer data 
were obtained using an Inconel heater sheet attached to 
the surface and coated with liquid crystals. Visual and 
quantitative flow field results using particle image velo- 
cimetry were also obtained for a plane at midchannel 
height for a Reynolds number of 45,000. The flow was 
seeded with polystyrene particles and illuminated by a 
laser light sheet. 

Computational results were determined for the 
same configuration and at matching Reynolds numbers; 
these surface heat transfer coefficients and flow velocities 
were computed with a commercially available code. The 
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experimental and computational results were compared. 
Although some general trends did agree, there were 
inconsistencies in the temperature patterns as well as in 
the numerical results. These inconsistencies strongly 
suggest the need for further computational studies on 
complicated geometries such as the one studied. 

Nomenclature 

A heated area, m 2 

h heat transfer coefficient, W/m 2 K 

Q heat rate, W 

r recovery factor 

T temperature, °C 

U w free-stream velocity, m/sec 
v kinematic viscosity, m 2 /sec 
p density, g/cm 3 

r w wall shear stress, N/m 2 

Subscripts: 

e electrical heat input 

L loss 

LC liquid crystal 
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R recovery 
s static 

t total 

Introduction 

A continuing objective in gas turbine technology is 
higher engine efficiency. One method of obtaining higher 
efficiency is to use higher engine operating temperatures. 
Such higher turbine inlet temperatures often require the 
use of internal cooling in the blades and vanes to main- 
tain surface temperatures thp.t the material can with- 
stand for sustained periods of time. When internal 
cooling is used, it is important that it be used in the 
most effective manner. Effective design of internal cool- 
ing passages requires an understanding of the flow and 
temperature patterns that occur, and it is greatly 
enhanced by the availability of dependable, validated 
computer codes. 

An experimental study was performed to obtain 
quantitative visual information on flow and heat trans- 
fer in a test section that had several significant features 
of the internal passage of an experimental cooled radial 
turbine blade described in Ref. 1. The objective of this 
study was to generate a set of experimental heat transfer 
and flow field data that could be used to validate com- 
puter codes for internal cooling systems. In recent years, 
there have been several attempts to numerically model 
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internal flow and heat transfer with computer codes. ’ 
Reference 2 points out that it is the combination of a 
complete experimental database and accurate, efficient 
analytical tools that will ultimately shorten design times 
and optimize future turbines. Reference 4 further em- 
phasizes the importance of coordination and cooperation 
between experimental and computational researchers in 
aeropropulsion systems. Reference 5 reports on a code 
prediction for a rotating square duct that was verified 
by experimental data. Through such coordination, the 
results of the present study will be used to validate 
internal cooling codes for more complicated geometries. 

Entrance velocity profiles and free-stream turbu- 
lence intensity were measured to establish upstream 
boundary conditions for the computations. Surface heat 
transfer coefficients were measured at entrance Reynolds 
numbers of 45,000, 335,000, and 726,000 (based on the 
entrance hydraulic diameter). The heat transfer data 
were obtained by using a uniform Inconel heater sheet 
attached to the test section floor and coated with ther- 
mochromic liquid crystals. Liquid crystals, by virtue of 
their property of changing color with temperature, pro- 
vide a means of measuring surface temperature and visu- 
alizing thermal patterns at any desired location. The 
local heat transfer coefficients were calculated by meas- 
uring the electrical input to the heater sheet under flow 


conditions and using a heat balance. Reference 6 de- 
scribes the method in more detail. 

Quantitative visual flow field measurements from 
particle image velocimetry are also presented for a plane 
at midchannel height and a test section entrance 
Reynolds number of 45,000. The flow field was seeded 
with l-/xm polystyrene particles and illuminated by a 
laser light sheet. This flow field work is an extension of 
work done in conjunction with the University of Akron 
where the procedure and the computerized data reduc- 
tion method were developed. Details of this method are 
given in Ref. 7. The present report includes a much 
larger extension of the flow field covered in Ref. 7 such 
that it can be compared with the heat transfer results. 

Computational results for the same configuration 
were obtained with a commercially available code. This 
code was used to calculate surface heat transfer coeffi- 
cients and flow field velocities at Reynolds numbers 
matching some of the experimental results. The experi- 
mental and computational results were compared. 
Although some general trends did agree, there were 
inconsistencies in the heat transfer patterns and magni- 
tudes. These inconsistencies demonstrate the need for 
further development of computational studies on compli- 
cated geometries such as the one investigated. 

Apparatus 

Figure 1 is a schematic of the complete test facil- 
ity. Room temperature air was drawn through a tunnel 
by a vacuum exhaust system. The two-piece test section 
(Fig. 2) was made of clear acrylic plastic and had a 
passage height of 5.08 cm. The entrance section was 
15.24-cm wide, and the main test section was divided by 
a partition into two passages 12.7-cm wide. One passage 
was open, and the other passage contained three rows of 
0.95-cm diameter pins in a staggered array. These pins 
were spaced 2.54-cm apart laterally in rows also spaced 
2.54-cm apart. This configuration captures the signifi- 
cant features of the cooled radial turbine blade shown in 
Fig. 3: 1 namely, the flow split created by a partition, 
the sharp corner, the cylindrical pins for turbulence gen- 
eration, and leading-edge stagnation regions. Entrance 
velocities were measured by the differential between a 
wall static pressure tap and a total pressure probe at the 
entrance. Pressure, temperature, and flow data were 
recorded by a commercial data-acquisition system. 
Photographic data were taken with one or more video 
cameras. 

For the heat transfer tests, a 0.025-mm thick 
Inconel sheet with a known uniform electrical resistance 
was attached to the floor of both the entrance and main 
test sections with a double-coated adhesive material. 
Each sheet had copper bus bars spot welded to it, was 
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coated with black lacquer, then spray painted with 
liquid crystals. A cross section of the test surface 
(Fig. 4) shows the 2.5-cm-thick foam insulation that 
prevented heat loss through the floor. For the heat 
transfer tests, the cylindrical pins were made of copper 
in order to simulate an actual blade where the surface 
heat is transferred to the pins. The pins were not 
attached to the surface but rested on the surface simply 
by gravity. 

The test section for the flow visualization tests was 
identical in configuration to the heat transfer test sec- 
tion but without the liquid-crystal heater sheets and 
without the foam insulation. Another difference is that 
the pins were made of clear acrylic plastic instead of 
copper. This allowed the laser light sheet to partially 
pass through the pins to illuminate the flow field. 

Procedure 

Liquid-Crystal Calibration 

The liquid-crystal layer was calibrated for color- 
temperature correspondence by a thermocouple attached 
to each heated Inconel sheet. Dark blue corresponded to 
the highest temperature; blue, green, yellow, and 
reddish-brown corresponded to decreasing levels of tem- 
perature. Yellow occurred over the narrowest tempera- 
ture band of 0.06 °C and was therefore the color used 
for temperature measurement. The yellow color location 
could be changed by varying the electrical heat input to 
the sheet under flow conditions. During calibration, 
when the yellow color was at the same location as the 
thermocouple, the temperature reading was 37.8 °C. 
This was, therefore, the calibration temperature. 

Experimental Measurements 

Entrance Velocity Profile and Turbulence Inten- 
sity . Velocity profiles were measured near the upstream 
end of the entrance test section by a surface static 
pressure tap and a traversable total pressure probe. As 
measurements were made, the probe was moved in small 
increments from the floor surface to a location just 
above the midchannel height position. It was physically 
impractical to measure the profiles at higher positions. 
Because the entrance passage was symmetrical about the 
midchannel height position, it was assumed that the 
velocity profile from midheight to the test section ceiling 
was identical to the measured profile from midheight to 
the floor. Free-stream turbulence intensity was measured 
by hot wire anemometry. All the measurements were 
made at entrance Reynolds numbers of approximately 
45,000, 335,000, and 726,000. 

Heat Transfer . Heat transfer tests were performed 
at the same entrance Reynolds numbers used for the 


other measurements: 45,000, 335,000, and 726,000. 
These corresponded to a range of velocities from 9.14 to 
155.5 m/sec. The operating procedure was initiated by 
bringing the tunnel to the velocity necessary to produce 
the desired Reynolds number. Electric power was then 
supplied to both entrance and main test section Inconel 
sheets at a matching heat flux, raising the temperature 
of the liquid crystal to the point at which color patterns 
appeared. The location of the yellow color lines was 
varied over the surface by changing the power input 
while the free-stream velocity remained the same. Be- 
tween data recordings, typically about 15 min was given 
to allow thermal equilibrium to occur. Data were taken 
by four charged-couple device (CCD) color video cam- 
eras mounted high above the test section. One camera 
viewed the entire entrance section, another camera 
viewed the entire main test section, and the final two 
cameras gave closeup views of the cylindrical pin area. 
All four cameras were activated simultaneously to 
achieve time correspondence, and the data were recorded 
on Super VHS video tape. A photograph of the test sec- 
tion under flow conditions in shown in Fig. 5. The color- 
temperature patterns on the surface are clearly visible. 

Flow Visualization . Flow visualization tests were 
performed at only one entrance Reynolds number, 
45,000. In these tests, the flow was seeded with 1-fim 
polystyrene spheres that followed the flow field. These 
particles were suspended in a heated water-ethyl alcohol 
mixture and sprayed into the flow stream with a spray 
nozzle. By the time the mixture reached the test section, 
much of the liquid had evaporated. The particles were 
illuminated by a 7-W argon-ion laser beam that was 
focused, through a series of mirrors and lenses, to a sheet 
of light approximately 2 mm thick. The result was a 
thin, illuminated, horizontal slice of the flow field at 
midchannel height. 

A single black and white CCD video camera re- 
corded the flow, showing the trajectories of the particles 
as streaklines on each video frame. A light chopper was 
placed in the laser beam to measure the velocity of the 
flow field. This chopper, which consisted of a rotating 
disk with holes spaced such that the light was periodi- 
cally blocked, could chop light at a maximum frequency 
of 3000 Hz. The instantaneous velocity of the particles 
was represented by the distance between interruptions in 
the streaklines on any one video frame as illustrated in 
Fig. 6. Similarly, the acceleration of the particles was 
represented by the change in the distance between inter- 
ruptions in the streaks. Because of the high magnifica- 
tion required to see the particle traces, only one small 
area could be viewed at a time. A computer-aided design 
(CAD) software package was used to combine the data 
from approximately 300 individual frames so that all of 
the information could be presented in one composite 
picture. 
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Data Reduction 


Heat Transfer . Surface heat transfer coefficients 
were calculated as follows: 

* _ Qe - Ql 
A(T lc - T r ) 

The local heat transfer coefficient h was calculated at 
the location of the calibrated color band (an isotherm— 
which in this case also represents a uniform heat transfer 
coefficient). The heat energy Q e supplied to the heater 
sheet was calculated from the measured voltage across 
the sheet and the current through a shunt resistor in ser- 
ies with the sheet. The heat loss Q L was the calculated 
radiation from the surface. The insulation under the test 
cell floor prevented back side heat loss, and other heat 
losses were considered negligible. The area A was the 
measured area of the Inconel heater sheet, and the tem- 
perature of the test surface T LC was the calibrated 
liquid-crystal temperature. The recovery temperature 
T r of the entrance free-stream air was calculated as 

Tr = T g + r(T t - T.) 

where T g and T t are the static and total air tempera- 
tures, respectively, and r is the recovery factor, defined 
as the cube root of the Prandtl number. 

An uncertainty analysis based on the method de- 
scribed in Refs. 8 and 9 was performed on the calcula- 
tion of experimental heat transfer coefficients. With this 
method, the maximum expected uncertainty values in 
the calculation of heat transfer coefficients were 14.0, 
3.7, and 3.8 percent for Reynolds numbers of 45,000, 
335,000, and 726,000, respectively. The relatively high 
percentage heat transfer uncertainty at the lowest 
Reynolds number resulted from a high uncertainty in 
flow velocity because of the lower accuracy of the 
pressure- measuring devices at very small pressure differ- 
entials. At all Reynolds numbers, errors were also possi- 
ble in reading the exact location of the isotherms. 

The heat transfer video data were reduced with a 
personal computer equipped with a commercial video 
frame grabbing board. A computer program was written 
to manually digitize (trace) the yellow color bands from 
a video image grabbed from the video tape, and the 
digitized image was stored in a computer. The digitized 
locations along with the measurements of surface heat 
flux, air temperature, and free-stream velocity allowed 
the calculation of heat transfer coefficients. 


Flow Visualization . The method for reducing the 
flow visualization data is described in detail in Ref. 7. In 
the present study, the region of the test section mapped 
with this technique was greatly increased to include an 
area from 8.9-cm upstream of the front row of pins to 
3.8-cm downstream of the last row of pins. 

The uncertainty involved in determining velocity 
is largely dependent on the accuracy of the operator dur- 
ing the digitizing process. Reference 7 reports that maxi- 
mum expected errors may vary from 4 to 8 percent. 

Computational Solutions 

Computational solutions of the flow field and heat 
transfer coefficients along the floor of the test section 
were obtained with the commercially available code, 
ST AR-CD. 10,11 This code solves the fully three- 
dimensional, turbulent compressible Navier-Stokes equa- 
tions on an unstructured mesh with local mesh refine- 
ment capabilities. The governing equations of mass, 
momentum, energy, turbulent kinetic energy, and dissi- 
pation Me discretized by the finite volume method, 
and a semi-implicit formulation is utilized (SIMPLE 
scheme). 12 

The computational mesh for this analysis is shown 
in Fig. 7. Because computer resources were limited, the 
computational grid was limited to about 250,000 cells, 
with finer mesh generated in the pin region. Boundary 
conditions were set at the entrance plane and at the exit 
of the computational domain to match those of the 
experimental values at Reynolds numbers of 45,000, 
335,000, and 726,000. Because of difficulties encoun- 
tered trying to input the measured velocity profile into 
the code, a constant velocity profile was used at the 
entrance. This profile was the free-stream velocity meas- 
ured at the entrance plane of the test section. In addi- 
tion, the average experimentally measured turbulence 
level of 2.3 percent was specified. A constant pressure 
was assigned at the exit, and the walls were specified to 
be at a constant temperature. 

Steady-state calculations with the upwind differ- 
ence scheme were chosen for this analysis. The k-e tur- 
bulence model was used, and heat transfer coefficients 
were obtained with the Reynolds analogy, which relates 
the convective heat transfer coefficient to the shear 
stress of the flow at the wall. A Unix-based, Silicon 
Graphics’ Iris workstation with four central processing 
units (cpu’s) in parallel was used to run the analysis and 
to pre- and post-process the data. The cpu time for 
obtaining a convergent solution ranged from 29 to 55 hr. 
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The converged solution of the flow provided veloci- 
ties, pressures, temperatures, density variations, and 
various other parameters throughout the channel, as 
well as heat transfer coefficients on the surface. 

Results and Discussion 
Experimental Results 

Entrance Velocity Profiles . Velocity profiles at the 
entrance are shown in Fig. 8 for the three entrance 
Reynolds numbers of 45,000, 335,000, and 726,000. 
Fully developed flow had not occurred, and as expected, 
there was a slight decrease in boundary-layer thickness 
as velocity was increased. The boundary-layer thickness, 
defined as the point where the velocity is 99 percent of 
the free-stream value, was 1.55 cm at a Reynolds num- 
ber of 45,000, 1.40 cm at 335,000, and 0.94 cm at 
726,000. 

The velocity data were also plotted logarithmically 
in terms of dimensionless coordinates (Fig. 9). This plot 
was made to show how the boundary layer compared 
with turbulent boundary layers normally encountered in 
turbine flow. Included in the figure is the calculated 
theoretical line for a fully turbulent boundary layer. The 
equations used were taken from Ref. 13. The experimen- 
tal data match the theoretical line fairly well. 

Free-Stream Turbulence Intensity . Free-stream 
turbulence intensity was measured with a hot wire ane- 
mometer at midchannel height. At a Reynolds number 
of 45,000 the turbulence intensity was measured as 

2.8 percent, at a Reynolds number of 335,000 it was 
1.2 percent, and at a Reynolds number of 726,000 it was 

2.9 percent. The average intensity was 2.3 percent. It 
has not been determined why the turbulence intensity 
was so much lower at the medium value of Reynolds 
number. 

Heat Transfer Measurements . Quantitative heat 
transfer coefficients on the floor of the test section were 
mapped for three different entrance Reynolds numbers. 
A typical example of the results is shown in Fig. 10. 
Parts (a) to (c) are actual video screen photographs of 
the surface at the lowest Reynolds number of 45,000 and 
at three of the various electrical energy inputs. The 
figure includes closeup views of the color-temperature 
patterns around the pins. Figure 11(a) shows the corre- 
sponding map of heat transfer coefficients calculated 
from these isotherms, and Fig. 11(b) is a closeup of the 
area around the pins. Similar heat transfer maps are 
shown in Figs. 12 and 13 for Reynolds numbers of 
335,000 and 726,000, respectively. 

Generally, the heat transfer patterns were similar 
for all Reynolds numbers. The heat transfer in the 


entrance region decreased going downstream as shown in 
Figs. 11(a), 12(a), and 13(a). Areas of low heat transfer 
are seen in the front left (relative to the flow stream) 
corner of the main test section where the flow is recircu- 
lating. The lowest values of heat transfer occurred in 
this region at all three flow conditions as shown in these 
same figures. In the unobstructed passage (no pins) there 
is a similarity in the pattern of the heat transfer at all 
three flow conditions with a high value heat transfer 
coefficient coming off the partition leading edge (proba- 
bly produced by the accelerating flow and high turbu- 
lence created in this area, as will be discussed later). The 
heat transfer coefficients decreased as flow progressed 
downstream. Close to the partition, the coefficients were 
even smaller, probably because of a flow separation in 
this area. From the closeup views of the pin region 
(Figs. 11(b), 12(b), and 13(b)), the cooling effectiveness 
of the pins is evident by the color patterns immediately 
downstream of the pins and the corresponding high val- 
ues of heat transfer coefficients. Further downstream of 
the pins, the heat transfer coefficients decreased. At 
some locations, the experimental contours were broad 
bands instead of thin lines, indicating very gradual sur- 
face temperature gradients. 

Flow Field Measurements . Flow field measure- 
ments were made at an entrance velocity of 9.14 m/sec, 
which corresponded to the lowest Reynolds number of 
45,000. Attempts were made to obtain results at the 
higher Reynolds numbers, but these attempts were un- 
successful with the equipment that was available. As the 
velocity of the flow field was increased, it was more dif- 
ficult to detect individual traces on the video screen. 
The maximum measurable particle velocity was further 
limited by the maximum frequency at which the light 
chopper could operate, which was 3000 Hz. As previ- 
ously shown (Fig. 6), the velocity was determined by the 
distance between interruptions in the streaks shown in 
any one video frame. At higher velocities, streaks would 
appear in a single frame with either one interruption or 
no interruptions at all, making it impossible to deter- 
mine the velocities or accelerations. 

Figure 14(a) is a video picture of the particles as 
they flowed between a pin and the partition leading 
edge. Figure 14(b) shows the data sifter being processed 
into velocity vectors. The numbers on the diagram iden- 
tify each particular trajectory, and the magnitude of the 
velocity is indicated by the length of the arrows. 

Figure 15 is a composite view of about 300 video 
frames, giving an overall view of the flow field in the 
area investigated. Velocity vector diagrams similar to 
the one shown in Fig. 14(b) can be generated anywhere 
within this area from the data collected in this study. 
Figure 15 reveals a few features that help explain some 
of the heat transfer results. A turbulent area appears 
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immediately downstream of the pins, which corresponds 
to the high heat transfer area seen in the maps. In the 
unobstructed passage where the flow comes off the parti- 
tion, there seems to be a region of turbulent flow bor- 
dered by an area of smooth accelerated flow. This 
corresponds to the high values of heat transfer seen on 
the heat transfer maps. In the front left corner of the 
main test section, the flow mapping appears to indicate 
a highly turbulent, high velocity area; but actual obser- 
vation of the flow showed this area to have very slow, 
three-dimensional flow, making it difficult to properly 
visualize the flow in a thin slice of light. The seeded air 
seemed to just circle slowly around in a lazy manner, 
which would lead to the low heat transfer and high sur- 
face temperatures indicated in the heat transfer tests. 

Computational Results 

Heat Transfer Computations . Computational solu- 
tions for the surface heat transfer coefficients were 
obtained for all three Reynolds numbers. The results are 
shown in Fig. 16 for a Reynolds number of 335,000 for 
comparison with the experimental results. Figure 16(a) 
shows a decreasing value of heat transfer coefficient in 
the entrance section. Low values were seen in the recir- 
culation region as expected. In the area downstream of 
the pins, the isotherms occurred in a generally longitudi- 
nal direction. Figure 16(b) shows a closeup view of the 
area around the pins. Low values were calculated in the 
wake of the pins. 

Flow Field Computations . Computational solutions 
for velocities at midchannel height of the flow field are 
shown in Fig. 17. The color chart designates the magni- 
tude of the velocities. These solutions were obtained for 
a Reynolds number of 45,000, which corresponds to an 
entrance velocity of 9.14 m/sec. As the flow approached 
the pins and the partition area, it accelerated into the 
unobstructed passage, whereas the flow going through 
the pins decelerated. The recirculation region shows very 
low velocities in the range of 2 to 4 m/sec. The flow in 
the pin area does not seem to be as turbulent as ex- 
pected. Evidently, the grid around the pins was not fine 
enough to show the complex mixing that occurs in the 
wake regions. 

Computational-Experimental Comparisons 

Heat Transfer Comparisons . There is agreement in 
some of the general trends between the experimental and 
computational results (for example, the decreasing heat 
transfer in the entrance region), but there is disagree- 
ment in the numerical values around and downstream of 
the pins. Because the Reynolds analogy was used in the 
analysis, this calculation would not be expected to be 
valid around the pins. In addition, the computational 
map, which shows longitudinal, fingerlike contours 


downstream of the pins, differs from the experimental 
maps, which showed almost constant heat transfer lines 
across the channel. Figure 16 shows high heat transfer 
coefficient values in the high acceleration region between 
the partition and the recirculation area, whereas Figs. 11 
to 14 show high heat transfer coefficient values in the 
high turbulent flow coming off the partition. The use of 
a finer mesh grid near the floor, where the heat transfer 
is calculated, might have produced more compatible 
results. In the recirculation region, the calculation shows 
a low heat transfer coefficient, which is confirmed by the 
experiment. 

Flow Field Comparisons . Because experimental 
flow field data were obtained at only one entrance flow 
velocity, 9.14 m/sec, a comparison between computa- 
tional and experimental results could be made only at 
the one corresponding Reynolds number of 45,000 and 
only at one small area at a time. The area chosen was 
the region near the leading edge of the partition and the 
nearest pin. A close-up view of the flow in this area is 
shown in Fig. 18. (This view is at midchannel height.) 
In the area upstream of the pins, there is good agree- 
ment with experimental velocities at about 7.5 m/sec. 
However, in the area between the partition and the pin, 
the code showed the flow decelerating to 5.5 m/sec, 
which is lower than experimental values seen in 
Fig. 14(b). The calculated values for the flow field show 
the general trends, but failed to predict the levels and 
details around the pins. The grid that was used to com- 
pute the flow field around the pins consisted of 16 cells 
around the pins and 20 cells deep — evenly spaced from 
top to the bottom of the channel. The number of the 
grid points was chosen arbitrarily to maintain a reason- 
able size grid and to obtain a solution in a timely 
manner. Using a denser grid or assuming symmetry 
about the midchannel height and analyzing only the 
lower half with the same number of grid points might 
have generated a better flow field calculation. Such 
refinements will be attempted at a later date. 

Concluding Remarks 

The experimental investigation revealed the tem- 
perature patterns that occur throughout a test section 
with some significant features of an internal passage of 
a cooled turbine blade. Areas of high and low heat 
transfer were identified, and quantitative heat transfer 
data were obtained over a range of flow conditions. The 
cooling effectiveness of the cylindrical pins was demon- 
strated. Areas of recirculating and separating flows and 
their effect on surface temperature patterns were evi- 
dent. The experimental flow visualization data generally 
supported the heat transfer data in terms of explaining 
some of the temperature patterns. Two benchmark sets 
of data were produced, heat transfer and flow field, 
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which are now available for internal cooling code 
validation. 

Computational solutions for the flow field and 
surface heat transfer were obtained with a commercially 
available code, STAR-CD. These solutions showed gen- 
eral trends of the flow field in the channel. The mag- 
nitude of the velocity field around the pins was lower 
than the experimental values. The code-predicted heat 
transfer results agreed with some of the general trends of 
the experimental data, but there were differences in the 
numerical values of the surface heat transfer coefficients 
on most of the channel floor. In the pin area, the calcu- 
lated heat transfer coefficients downstream of the pins 
showed trends opposite to those of the experiment. The 
grid size of 250,000 cells was evidently insufficient to 
resolve the flow field in this area as well as the heat 
transfer. 

STAR-CD, which is typical of currently available, 
general purpose, computational fluid dynamics (CFD) 
codes, showed the limitations that exist at the present 
time in predicting the flow field and heat transfer in 
complex geometries. The difficulties encountered in 
attempting to match computational simulation and ex- 
perimental data demonstrate the need for further devel- 
opment of CFD codes. 
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Figure 1 . — Branched duct tunnel. Note: Only one camera used for flow-visualization tests. 
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Figure 3. — Cooled radial turbine blade (Ref. 1). 
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Figure 5— Test section underflow conditions; Reynolds number, 45,000. 
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Figure 6. — Instantaneous velocity from 1 /30-sec video frame 
using 3000-Hz light chopper. 
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Figure 8.— Entrance velocity profiles. 
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Figure 9.— Dimensionless boundary-layer profile. 
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(a) Low-heat input; main section. 

Figure 10. — Surface color-temperature patterns; Reynolds number 45,000. 
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(b) Medium heat input. 
Figure 1 0. — Continued. 
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(c) High heat input. 

Figure 1 0. — Concluded. 
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(a) Entire test surface. 



(b) Close-up of pin area. 


Figure 1 1 .—Experimental heat transfer coefficients; Reynolds number, 45,000. (Coefficients given in W/m 2 *K.) 
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Separation region 




(b) Close-up of pin area. 


Figure 1 2.— Experimental heat transfer coefficients; Reynolds number, 335,000. (Coefficients given in W/m 2 «K.) 
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(a) Entire test surface. 



(b) Close-up of pin area. 


Figure 13. — Experimental heat transfer coefficients; Reynolds number, 726,000. (Coefficients given in W/m 2 *K.) 
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(a) Video picture. 



Figure 1 4. — Particle trajectories; Reynolds number, 45,000. 
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Figure 1 5. — Composite of digitized video data; Reynolds number, 45,000. (300 video frames combined.) 
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Figure 17. — Computational flow field map; Reynolds number, 45,000. 
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